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Abstract 

The purpose of this paper is to present a new mathematical model for 
the deformation of thin Cosserat elastic plates. Our approach, which is 
based on a generalization of the classical Reissner plate theory, takes into 
account the transverse variation of microrotation of the plates. The model 
assumes polynomial approximations over the plate thickness of asymmet- 
ric stress, couple stress, displacement, and microrotation, which are con- 
sistent with the elastic equilibrium, boundary conditions and the consti- 
tutive relationships. Based on the generalized Hellinger-Prange -Reissner 
variational principle and strain-displacement relation we obtain the com- 
plete theory of Cosserat plate. We also proved the solution uniqueness for 
the plate boundary value problem. 
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1 Introduction 

The well known classical bending theory of elastic plates [5], [7], [T7], was first 
presented by Kirchhoff in his thesis (1850) and is described by a bi-harmonic 
differential equation i2],[T3. The usual assumption of this theory is that the 
normal to the middle plane remains normal during deformation. Thus the theory 
neglects transverse shear strain effects. A system of equations, which takes into 
account the transverse shear deformation, has been developed by E. Reissner 
(1945) [13], [g. 

One of the advantages of Reissner's model is that it is able to determine 
the reactions along the edges of a simply supported rectangular plate, where 
classical theory leads to a concentrated reaction at the corners of the plate. 
The Reissner theory has been applied to thin walled structures with moderate 
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thickness. The study of the relationships between these two models has proved 
that the solution of the clamped Reissner plate approaches the solution of the 
Kirchhoff plate as the thickness approaches zero [T] and that the maximum 
bending can reach up to 20% for moderate plate thickness [5]. The numerical 
calculations of bending behavior of the plate of moderate thickness, (TH] show 
high level agreement between 3D and Reissner models. More remarks on the 
history of the modeling of classic linear elastic plates can be found in [6] , [16] , 

In order to describe deformation of clastic plates with microstructure that 
possess grains, particles, fibers, and cellular structures [TD], [Tl]. A. C. Eringen 
(1967) was the first to propose a theory of plates in the framework of Cosserat 
(micropolar) Elasticity [3] . His theory is based on a direct technique of integra- 
tion of the Cosserat Elasticity. The Eringen plate theory does not consider a 
transverse variation of the microrotation over the thickness, which might be nec- 
essary for rather thick plates under vertical load and pure twisting momentum. 
In order to develop a theory of plates, which can be used for thin wall structures 
with moderate thickness, we propose to use the classic Reissner plate theory as 
a foundation for the modeling of Cosserat elastic plates. Our approach, in addi- 
tion to the traditional model, takes into account the second order approximation 
of couple stresses and the variation of three components of microrotation in the 
thickness direction. 

2 Micropolar (Cosserat) Linear Elasticity 
2.1 Fundamental Equations 

Before proceeding some notation convention should be explained. We use the 
usual summation conventions and all expressions that contain Latin letters as 
subindices are understood to take values in the set {1, 2, 3}. When Greek letters 
appear as subindices then it will be assumed that they can take the values 1 or 
2. 

The Cosserat elasticity equilibrium equations without body forces represent 
the balance of linear and angular momentums of micropolar elasticity and have 
the following form [3]: 

divcr= , (1) 

£ • cr +div/x = 0, (2) 

where the quantity a = {<Jji} is the stress tensor, fi = {/iji} the couple stress 
tensor, e = {sijk} is the Levi-Civita tensor, where Sijk equals 1 or -1 accord- 
ing as {i,j,k) is an even or odd permutation of 1,2,3 and zero otherwise, and 
£ ■ cr = {SijkfJjk} ■ 

The constitutive equations can be written in the form [12j : 

= (M + Mc)7 + (M-Mc)7^ + A(tr7)l, (3) 
M = (7 + e)x + (7-e)x^ + /3(trx)l, (4) 
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and the strain-displacement and torsion-rotation relations 

7 = (Vu)^ +£ • y> and x — (5) 

where quantities 7 and x i the micropolar strain and torsion tensors, u and 
(f the displacement and rotation vectors respectively, 1 the identity tensor, and 
/Lt, A are the symmetric, /i^,/9,7,e the asymmetric Cosserat elasticity constants. 
In the reversible form: 

7 = i^i' + ^i',)(T + ifi' -fi'^)a^ + X'itra)l, (6) 
X = (7' + e')/^ + (7'-e')M'^+/3(tr/x)l. (7) 

where ^' = ^, = ^, 7' = ^, ^ ± X' = ^^^^ and /3' = 

We consider a Cosserat elastic body Bq. In this case the equilibrium equa- 
tions (H]) - (H]) with constitutive formulas ([3]) - (g]) and kinematics formulas ([5]) 
should be accompanied by the following mixed boundary conditions 

u = Uo, ip = iPo on = dBa\dB„, (8) 
<T„ = (7 ■ n = (To, = fJ- ■ ^ = f^o on 02 = dB^, (9) 

where Uq, (p^ are prescribed on Qi, do and /x on ^2, and n denotes the outward 
unit normal vector to DBq. 

2.2 Cosserat Elastic Energy 

The strain stored energy Uc oi the body -Bo is defined by the integral [12]: 



Bo 



Uc= I W{j,x}dv, (10) 

where 



W{'y,x} = Y^ltjlji+ 2^kklnn (11) 

J + e 7 - e P 



then the constitutive relations ([3]) - (|4]) can be written in the form: 

a = C„[W\=\I^W&nd iJi^C^[W\ = V^W. (12) 
The function W is positive if and only if JJj 

^l > 0, 3A + 2/i>0, 

7 > 0, 3/3 + 27>0, (13) 
/i^ > 0, -I- /i^ > 0, e > 0. 



3 



The following conditions [3] for the Cosserat elastic energy 



/i > 0, 3X + 2fi> 0, 

7 > 0, 3/? + 27>0, (14) 
Mc > 0, e > 0. 



are enough to provide the uniqueness of static problems. 
For future convenience, we present the stress energy 



Uk = ${cr,/x}du, 



where 



The reversible constitutive relation ([6]) - ([7]) can be also written in form: 

7 = K,H = ^, X = KJ/x] = ^. (16) 

The total internal work done by the stresses <t and /x over the strains -y and 
X for the body So [E] is 



f7 = / [cr • 7 + /X • x] 



(17) 



and 

U = C/k= Uc 
provided the constitutive relations ([2]) - Q hold. 

2.3 The Generalized Hellinger-Prange -Reissner (HPR) 
Principle 

The HPR principle [S] in the case of Cosserat elasticity states, that for any set A 
of all admissible states s ~ [u, if, 7, x, cr, /x] that satisfy the strain-displacement 
and torsion-rotation relations ([5]), the zero variation 

(56 (s) = 

of the functional 

9(5) = Uk - I [(T--f + ii-x]dv (18) 

JBo 

+ [(T„ • (u - Uo) + (¥> - (Po)] da+ [a-o-u + nio- <fi] da 

J Qi J g2 
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at s £.4 is equivalent of s to be a solution of the system of equilibrium equations 

([1]) - constitutive relations (HI) - O, which satisfies the mixed boundary 
conditions ([5]) - ©■ The proof is similar to the proof for HPR principle for 
classic linear elasticity 



3 The Cosserat Plate Assumptions 

In this section we formulate our stress, couple stress and kinematic assumptions 
of the Cosserat plate. The set of points P = {T x [-h/2, h/2]} UTUB forms 
the entire surface of the plate and {r„ x [—h/2, h/2]} is the lateral part of the 
boundary where displacements and microrotations are prescribed. The notation 

= r\r„ of the remainder we use to describe the lateral part of the boundary 
edge {To- x [—h/2, h/2]} where stress and couple stress are prescribed. We also 
use notation Pq for the middle plane internal domain of the plate. 

In our case we consider the vertical load and pure twisting momentum 
boundary conditions at the top and bottom of the plate, which can be writ- 
ten in the form: 



cr33{xi,X2,h/2) = a\xi,X2), a33{xi,X2,-h/2) 

a3f3{xi,X2,±h/2) = 0, 

M33(2;i,a;2,ft./2) = ^i*-{xi,X2), ^i33{xi,X2, -h/2) 

M3/3(a;i,a;2,±ft-/2) = 0, 

where (xi,X2) £ ^o- 

3.1 Stress and Couple Stress Assumptions 

Our approach, which is in the spirit of the Reissner's theory of plates [13], 
assumes that the variation of stress (Jki and couple stress /.t^.; components across 
the thickness can be represented by means of polynomials of X3 in such a way 
that it will be consistent with the equilibrium equations ([T]) and ([2]). First, as 
it is assumed in the standard theory of plates, we use expressions for the stress 
components in the following form [18] : 

= nap{xi,X2) + ^Cz'TTiafiixi, X2) , (23) 

where ^3 = j^x^, and a,/3 £ {1,2}. The only difference between our assump- 
tions and those of Reissner' [T3] is that the functions Uap and nriap are not 
symmetric. Based on (j23p and by means of the first two equations of written 
in the component form stress equilibrium ([T]) 

we obtain for the shear stress components 



^a\x^,X2),m 
(20) 

= ^'(^1, 2:2), (21) 
(22) 
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cr3f3 ^ qi3{xi,X2) (l - Cl) , (24) 

It is natural to assume that the expressions for the remaining shear stress 
component are in the form similar to ([M)) . i.e. 

<^l33^qfi{xi,X2){l-Cl) ■ (25) 

Here, as is usual for the asymmetric elasticity, the functions qp^q*^. can be 
different. 

Substituting equations (|25p in the remaining equilibrium differential equa- 
tion for stress 

f^jsj = 

we obtain the expression for the transverse normal stress 

(^33 = C3 (^C3 - l) k*{x^,X2) + m*(xi,X2). (26) 

The next step is to accommodate approximations (|26p to the boundary con- 
ditions (|19p . By direct substitution to ^9]) it easy to obtain that 

3 /1. 3 



^33 = -^ (^3C3-C3jP + ^0, (27) 

where p = (t*(xi,X2) — (T''(a;i, X2) and ctq = 5 (c*(j:i, 0:2) + cr''(a;i, 0:2)) satisfy 
the boundary condition requirements. We note that expression ([27| is identical 
to the expression of 0-33 given in [13] in the case of cr'' = 0. 

It is also assumed that the couple stress /i^^ should have expression similar 
to the shear stress a^p expressions (|24|) - ([25l) : 



Ma/3 = (1 - C3) rci3{xi,X2). (28) 
Finally we assume that couple stress /i^3 expression is similar to aap (|23p : 

^^P3 = C7iSl{xi,X2) +rnlj{xi,X2). (29) 
Note that the first two equations of ([2]) can be written in the form 

e/jj-feCTifc + lJ.ji3,j = 0, (30) 

and substituting the couple stress ([^5]) in (|30p and taking into account ([M]) and 
([25]) we obtain the expression for the transverse shear couple stress: 



M3/3 = - C3^ Si3{xi,X2) + TO^(xi, 0:2). (31) 

Substituting ((3T|) to boundary conditions ((2T|) we obtain that 

si3{xi,X2) = and m/3(xi,a;2) = 0, 
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i.e. the transverse shear couple stress 



= 0. (32) 

Now, substituting the couple stress ((29)) and stress ([23]) in the remaining differ- 
ential equation of the equilibrium of angular momentum ((2]) 

esjfecrjfe + /^j3,j = 0, (33) 
we obtain the transverse normal couple stress to be in the form: 

M33 = ^C3a*(2;i, 2:2) + C3^*(a;i, a;2) + c*{xi,X2). (34) 

The next step is to accommodate boundary conditions (PT|) to ([M)) . At this 
stage we restrict the form of which could allow us to determine couple 
stress /i33 directly from boundary conditions (|2ip . To this end we make an 
additional assumption that /igg must be a first order polynomial 

/^33 = C3^*(2;i,a;2) + c*(a;i,a;2). (35) 

This assumption is also consistent with the equilibrium equation (|33p and al- 
lows us to proceed as we did for the determination of transverse loading stress 
(j27p from the stress boundary conditions. Now boundary conditions (|2ip are 
sufficient to determine /X33, which must be of the form 

M33 = C3V + t, (36) 

where ?;(a;i, 2:2) = ^ (^*(xi, 2:2) — ^''(xi, X2)) and t{xi,X2) — \ (/i*(a;i, 0:2) -I- ^J^[xi 



3.2 Kinematic Assumptions 

The choice of kinematic assumptions is based on simplicity and their compati- 
bility with the constitutive relationships of stress and couple stress assumptions 
([3]). As in the standard theory of thin plates, it is assumed that Ua displace- 
ments are distributed linearly over the thickness of the plate 3 and that 1*3 
does not vary over the thickness of the plate, i.e. 



Ua = Ua{xi,X2) -^C,3Va{xi,X2), (37) 
"3 = w{xi,X2), 

The terms Va{xi^X2) in ()37p represent the rotations in middle plane. 

In order to accommodate the transverse microrotations to the constitutive 
relations ^ we propose the variation of microrotation with respect to X3 by 
means of the second and third order polynomials: 

= e°(a:i,a;2)(l-C3), (38) 
^3 = e°(a;i,:r2) + C3(l- JC3)03(a;i,X2). (39) 
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The constitutive formulas ([2]) - ([U motivate us to chose the forms ([55]) and ([5^ , 
which produce expressions for Lp^ p and (^3 3 similar to what we have for couple 
stress approximations pS)) . 

The functions 0° in and (|39p describe microrotation components in the 
middle plane of the plate and 03(2:1, X2) the slope at the middle plane. Thus, in 
the assumptions ([55)1 and the transverse variation effect of microrotations 
is not neglected. 

4 Specification of HPR Variational Principle for 
the Cosserat Plate 

The HPR variational principle for a Cosserat plate is most appropriately ex- 
pressed in terms of corresponding integrands calculated across the whole thick- 
ness. We also introduce the weighted characteristics of displacements, micro- 
rotations, strains and stresses of the plate, which will be used to produce the 
explicit forms of these integrands. 

4.1 The Cosserat plate stress energy density 

We define the plate stress energy density by the formula; 

<^{S)^^ j\{a,t,]d(:^. (40) 

Taking into account the stress and couple stress assumptions ([23l) - ()36|) and 
by the integrating $ {a, /x} with respect C3 in [—1,1] we obtain the explicit 
plate stress energy density expression in the form: 
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$(5) 



2ft/i(3A + 
A 

2/i/i(3A + 



12 



3 



(1 - ( + ^MlA + ^ (g„Q„ + 



5 10 



3A 



3 



5/i/i(3A + 2fi) 
1 1 



^^7MA + Ai)_^2 



A 



280Ai(3A + 2/i) ^ 2/x(3A + 2/x) 



/j^A+^i) 2 



7 + e 



2/i(3A + 2/i) ° /i7e [8 



Mx + ^s:s: + -{i-6^,)R 



20 



^7 



27(3/3 + 27) 
where the Cosserat stress set 



,2 ^ J£±7L,2 



where 



27(3/3 + 27)^ ' 67(3/3 + 27) 
5 = [M„^, Qo,, Ql^, Raf}, Sp, Naf}, M*] , 



(41) 



(42) 



Qo 



Rafi — 



<3„ — 



N, 



a/3 



2 .1 ^3 
f3a"C3 = -;r9c«) g, 



/I 2/1 

2 y ^Ma/jaCs = Y»^a/3, 

(^) /_/3M„3rfC3 = 



0-a3aC3 = ^9, 



dCs = hm* 



(43) 



Here Mn and M22 are the bending moments, M12 and M21 the twisting 
moments, Qa the shear forces, g* the transverse shear forces, i?ii and i?22 the 
micropolar bending moments, -R12 and -R21 the micropolar twisting moments, 
5* the micropolar couple moments, all defined per unit length, iVn and N22 are 
the bending forces, iVi2 and A^2i the twisting forces, M* the micropolar shear 
couple-stress resultants. 
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Then the stress energy of the plate P 




$(5)da, 



(44) 



where Pq is the internal domain of the middle plane of the plate P. 

4.2 The density of the work done over the Cosserat plate 
boundary 

In the following consideration we also assume that the proposed stress, couple 
stress, and kinematic assumptions are valid for the lateral boundary of the plate 
P as well. 

We evaluate the density of the work over the boundary r„ x [—h/2,h/2] 



Taking into account the stress and couple stress assumptions (|23|) - (|36|) and 
kinematic assumptions ([37]) - ([39|) we are able to represent Wi by the following 
expression: 



Wi - Sn-U -A/^*^ + Q*W + Rc^fll + 3*^3 + N^Ua, + M*nl, (46) 




(45) 



where the sets 5„ and U are defined as 



5, 



71 



and 



S* = 



In the above np is the outward unit normal vector to Fu, and 
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1 

Here are the rotations of the middle plane around W the vertical 

deflection of the middle plate, fl^ the microrotations in the middle plate around 
Xk axis, Ua is the in-plane displacements of the middle plane along Xa axis, fis 
the rate of change of the microrotation cp^ along X3. 

We also obtain the correspondence between the weighted displacement and 
the microrotations (|T7)) and the kinematic variables by applying ([57]) and (|39p 
in integration of expressions (j47p : 



*a = Vc{xi,X2), W ^W{X1,X2), (48) 
^1 - kiel{xi,X2), n3='^Q3{xi,X2), 

Ua = C/a(a;i,a;2), fi" = 63(2;i,a;2), 

where coefficients fci and ^2 depend on the variation of microrotations. Under 
the conditions ([55]) we have that fci = | and k2 — ^■ 

The density of the work over the boundary Fo- x [—h/2, h/2] 

h 

W2 = - y {cToaUa + rUoafa) ^0(^(3 

can be presented in the form 

W2 = So-U =Iioa^a + lio^W + MoMl + ACa^is + So,«[/« + Tosnl 
where 

Mapnp = Uoa, Rapnis = Moa, 

Q>a = Uo3, SlUa = ACg. (49) 
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Naf3n0 = E„, (50) 
M>„ = To3. (51) 



Now 71/3 is the outward unit normal vector to T^, and 
h h 



-1 



^o,a = 2/ <^o"^C3: Tf"o3 = ( 77 ) / C3 if^o3 - CsV) dCs- (52) 



We are able to evaluate the work done at the top and bottom of the Cosserat 
plate by using boundary conditions (fTU)) and (^11) 

(o'o3^t3 + »TT-o3<Po3) "srfa = {pW + vn°)da. 



TUB 

Pa 



4.3 The Cosserat plate internal work density 

Here we define the density of the work done by the stress and couple stress over 
the Cosserat strain field: 



h 

W3 = 2 / ^('^■7 + M-x)rfC3- 



(53) 



Substituting stress and couple stress assumptions - ([55]) and integrating 
expression ([53]) we obtain the following expression: 



(54) 

where £ is the Cosserat plate strain set of the the weighted averages of strain 
and torsion tensors 
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Here the components of £ are 

ea/3 = 



IJ^Cslc^pdC^, (55) 
f^lc.A^-C^)d(:s, (56) 



3 

3 

4 



1 



73c.(l-C')dC3, (57) 



3 

T3a = J^J ^CsXsadCs, (58) 



rc./3 = Ij\af3{^-C')d(:„ (59) 



1 



= ^lafidCs^ (60) 

Xaa^Ca- (61) 



1 



3q 



2 7-1- 

The components of Cosserat plate stram ([55|) - (|5l]) can also be represented 
in terms of the components of set U by the foUowing formulas: 





= 


+ £30/3^3, 


UJa 


= + 


£3af3^^, 






- £30/3^^/3, 


T 3a 






_o 


- QO 

— "/3,a; 






= U/3^a ■ 


■f £3a/3f^3j 




- fiO 





(62) 



We call the relation (p^ the Cosserat plate strain-displacement relation. 

5 Cosserat Plate HPR Principle 

It is natural now to reformulate HPR variational principle for the Cosserat plate 
P. Let A denote the set of all admissible states that satisfy the Cosserat plate 
strain-displacement relation and let 8 be a HPR functional on A defined 

by 

Q{s) ^Ulc- j {S ■ E-pW + vVtl)da + j Sn {U - Uo) ds + j So-Uds, (63) 
for every s = [U,£,S] G ^.Then 

Se{s) = 
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is equivalent to the following plate bending (A) and twisting (B) mixed prob- 
lems. 

A. The bending equilibrium system of equations: 



Mc0,a-Ql3 = 0, (64) 

Q*a,c.+P = 0, (65) 

Ra0,a + £.3/37 {Q^ " Q7) = 0' (66) 

S:,c. + e3PfM(3^ = 0, (67) 
with the resultant traction boundary conditions : 

Mafjnp = Uoa, Rapnp = Moa, (68) 

Q*na = IIoS, S^Ua = To3, (69) 

at the part and the resultant displacement boundary conditions 

*a = *oa, W = Wo, fll = 0°^, O3 = ^o3, (70) 

at the part r„. 

The constitutive formulas: 

. _ 3* _ + ^M.,- (71) 



II 6A 3A , . 

3(Mc + m) „^ , 3(/i,-M) ,^ ,^ 



(73) 



I 3/3 /3 

'^"'^^ ' 5/17(3/3 + 27) 27(3/3 + 27) ^' 

9$ 3(6 - 7) „ , 3(7 + 6) „ 

^ O-R/Ja 10/176 10/176 

a* 3(7 + 6) 
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B. The twisting equilibrium system of equations: 



^a/3,a - 0, (76) 
+ £3/57^/37 + « = 0, (77) 

with the resultant traction boundary conditions at To-: 

Napnp = Eq, (78) 
M>„ = m:„ (79) 

and the resultant displacement boundary conditions at T^- 

t/a = C/oa, = (80) 

The constitutive formulas: 



dNaa hfl{3X + 2fl)' 



2V(3A + 2Ai) 2^(3A + 2/i) 

a$ 7 + e 



We also represent the above constitutive relation in the compact form: 

where we call K, the compliance Cosserat plate tensor. 
Proof of the principle. The variation of 0(s) 

Seis) = J {{IC[S]-£)-SS-SS£+pSW + v6nl}da 

Po 

+ [ {SSr,-{U-Uo)+Sn-5U}ds+ [ So-SUds. 

We apply Green's theorem and integration by parts for S and SU to the 
expression: 

Js-6£da = J So6-U ds- J{iM^f3,c.-Qp)S^p + Ql^JW 

Po dPo Po 

+ {Rap.a + £3/37 (Q7 — Q7) Rap.a) 5fl^ 
+ {S*c,a. + £3/37^^/37) + N^p.aSUp 

+ {M'^^c + £3/37^/37) 5nl}da. 
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Then based on the fact that 5U and 5S satisfy the Cosserat plate strain- 
displacement relation (|62p . we obtain 



6@{s) = j {{]C[S] - 8) ■ 5S - S5£} da 

Po 

+ J {(M„0,„ - Qp) d^p + (Q:_„ + p) 6W 



'0 



+ {Ral3,a + £3/37 (Q7 — Q7) Rafl,a) Sflp 

+ {S*o,,a + £3/37^^/37) -^^^3 + N^p.aSUp + (M* „ + esp-^Np-^ + v) 5n°}da 



■SUds. 



iT„ Jr. 
If s is a solution of the mixed problem, then 

6e{s) =0. 

On the other hand, some extensions of the fundamental lemma of calculus 
of variations ,5] together with the fact that 14 and £ satisfy the Cosserat plate 
strain-displacement relation (|62p imply that 5 is a solution of the A and B 
mixed problems. 

Remark. In the case of T3„ = ^3,0 = we obtain that 5* = 0. and M is 
symmetric, i.e. Ma/3 — Mp^ and the corresponding constitutive relation is 



a$ _ 6 

The bending system for this case (j64l) -(|67 |) is reduced to the following: 

M^f3,a-Qp = 0, (84) 
Q:,a+P = 0, (85) 

RaP,a + £3^7 (Q7 — Qj) — 0. (86) 

We also notice 

and the case = is consistent with the requirements 

Q* = Qa and = W^a 
Thus we obtain the equilibrium system in the decoupling form: 

Ma(3.,c,-Q(3 = 0, (87) 

Qa,a+P - 0, (88) 
i?a/3,a = 0. (89) 
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6 Solution Uniqueness 



Here we prove that if there is a solution for the deformation of a Cosserat 
elastic plate, which satisfies the equilibrium equations (|64|) - ([67|) . (I76p - ((77)) . 
constitutive (|7ip - ([75]) , (|8ip - ([55]) and kinematics formulas (15^ with boundary 
conditions at To- and dTD]) and dM]) at then this elastic 

solution must be unique. We also assume that all functions and the plate middle 
plane region Pq satisfy Green - Gauss theorem requirements. 

The proof will be based on contradiction. Let us assume that the solution 
of the Cosserat plate is not unique in terms of the stresses and strains, i.e. there 
would be two different solutions of ([51)) - ([57)) and ([75)) - ([77)) . both of which 
satisfy the same boundary conditions ([^0)) and ([^^ at and ([70)1 and ([50)) at 
Fij. Due to linearity of the proposed model, the difference between these two 
different solutions is also a solution of the same system of equations with the 
following zero boundary conditions: 

M^^np = 0, Ro^pnp = 0, (90) 
Q*ana = 0, S*na = 0, 

Naf^n^ = , M*n^ = , (91) 

or 

■9^^o,w = 0, V p = o, n° = 0, ^3 = 0, u^ = 0. (92) 

It can be shown that for zero loads, the internal work U can be expressed 
by applying integration by parts as follows: 



U = J S-£da = J {{M^f^^^ + Q}W + Rc,pn'i + S*pn3 + N^f3Uc, + Mp°)^^ 

Po Po 

- {M^l3,a - Ql3) */3 - Ql. JV 

- {RaP,a + £3/37 (Q^ — Qf) Ra0,a) (93) 

- {Ml^ + £3/37^^/37) ^l}d'^- 

Taking into account Green's theorem, the equilibrium equations ()64p - ()67p and 
([76)) - ([77)) . expression ([93)) is reduced to the following line integral: 

U = <j> Sn-Uds = j, (Af„*„ + Q*W + Ro0l + S*^3 + NJJ^ + A?*r2°) ds = 0, 

(94) 

which vanishes because of the zero boundary conditions ([50)) - ([M)) . 

Using the constitutive equation ([7T)) - ([75)) . ()8ip - ()83p in a reversible form, 
the positive definite quadratic form strain energy density ()4ip can be represented 
in terms of the Cosserat plate strain set which components in this case 
should be zeros. Then from the Cosserat plate strain-displacement relation ()62p 
we obtain the system: 
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C/i,i = 0, U2,2 = 0, C/2,1 - 0° = 0, Ui,2 + fl°s = 0, 

+ = 0, w,2 -^1 = 0, = 0, *2,2 = 0, n^.i = 0, 03,2 = 0, 

*2,1 - fls = 0, *i,2 + fis = 0, *i - 0° 0, *2 + = 0, 

= 0, 0^,2 = 0, i^;,2 = 0, = 0, = 0, f]°,2 = 0, 

which has the foUowing solution: 

Ua = £30a^lx0 + U°, W = Saa^nlxff + W° , = £3a/3^^^, ^^3 = 0, 

where constant parameters C/", W'^ are the rigid translations of the middle 
plane, the component of the rigid rotation of the plane around and 
flz the slope of the rigid rotation around 0:3. Thus the difference between any 
two deformations and microrotations of the plate, having the same boundary 
conditions, represents changes of the plate as a rigid body. 

7 Conclusion 

We generalized Hellinger-Prange -Reissner (HPR) principle in order to derive 
the new equilibrium equations in the middle plane and constitutive relationships 
for the plate. The polynomial approximations of the variation of couple stress 
and micropolar rotations in the thickness direction in order higher than one 
allowed us, based on the generalized HPR principle, to project Cosserat 3D 
Elasticity equilibrium equations into the new form of equilibrium equations and 
the constitutive relations in the middle plane of the plate. 

8 Acknowledgement 

I express my gratitude to Krzysztof Rozga and Darrell Hajek for providing 
assistance in the presentation of the paper. 

References 

[1] Bathe K. J., Brezzi F.: On the convergence of a four node plate bending 

element based on Mindlin- Reissner plate theory and a mixed interpolation. 
The Mathematics of Finite Elements and Applications V ( Uxbridge,1984), 
Academic Press, London, 491-503, (1985) 

[2] Donnell L. H.: Beams Plates and Shells, McGraw-Hill, New York (1976) 

[3] Eringen A. C: Theory of micropolar plates. Journal of Applied Mathemat- 
ics and Physics, Vol. 18, 12-31, (1967) 

[4] Gauthier R.D., Jahsman W.E.: A quest for micropolar elastic constants. 
Journal of Applied Mechanics, 42, 369-374 (1975) 



18 



[5] Giirtin M. E.: The Linear Theory of Elasticity in Handbuch der Physik, 
Vol. VIa/2; C. Truesdell (editor), Springer- Verlag, 1-296 (1972) 

[6] Love A. E. H.: A Treatise on the Mathematical Theory of Elasticity, Dover, 
New York (1986) 

[7] Naghdi P. M.: The Theory of Shells and Plates, in Handbuch der Physik, 
Vol. VIa/2; C. Truesdell (editor). Springer- Verlag, 425-640 (1972) 

[8] NefF P.: A geometrically exact Cosserat-shell model including size ef- 
fects, avoiding degeneracy in the thin shell limit. Part 1: Formal dimensional 
reduction for elastic plates and existence of minimizers for positive Cosserat 
couple modulus. Cont. Mech. Thermodynamics, 16: 577-628 (2004) 

[9] Neff P.: The Cosserat couple modulus for continuous solids is zero 
viz the lineaqrized Couchy-stress tensor is symmetric. Preprint 2409,, 
httpI77www3.mathematik.tudarmstadt.de/fb/mathe/bibliothek/preprints.html 
(2005) 

[10] Neff P.: A finite elastic -plastic Cosserat theory for polycrystalls with grain 
rotations. Int J. Eng. Sci. 44574-594 (2006) 

[11] P. Neff and Forest S.: A geometricall exact micropolar model; for elastic 
metallic foams acounting for afhne microstructure. Modeling, existing of 
minimizers, identification of moduli and computational results. J. Elastic- 
ity, 87, 236-279, (2007) 

[12] Nowacki W.: Theory of Asymmetric Elasticity, Pergamon Press, Oxford, 
New York (1986) 

[13] Reissner E.: The effect of transverse shear deformation on the bending of 
elastic plates. Journal of Applied Mechanics, June, 69-77 (1945) 

[14] Reissner E.: On the theory of Elastic Plates, Journal of Mathematics and 
Physics, 23, 184-191, (1944) 

[15] Reissner E.: Reflections on the theory of elastic plates, Applied Mechanics 
Reviews, 38, 1453-1464 (1985) 

[16] Rossle A., Bischoff M., Wendland W., Ramm E.: On the mathematical 
foundation of the (l,l,2)-plate model. International Journal of Solids and 
Structures, 36, 2143-2168 (1999) 

[17] Timoshenko S. and Woinowsky-Krieger S.: Theory of Plates and Shells, 
McGraw-HiU (1959) 

[18] Wan F.Y.M.: Lectures Notes on Problems in Elasticty:II Linear Plate The- 
ory. Tech. Rep. No. 83-15, Institute of Apphed Mathematics, Univesity of 
British Columbia (1983) 



19 



